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ABSTRACT 


A barotropic and two- and three-layer baroclinic 
primitive equation prediction models are developed from a 
ten-layer model previously designed for short range fore- 
Casts in the tropics. Twenty-four hour predictions from 
the models are compared to evaluate any degradation in the 
forecasts as the vertical resolution of the models decreases. 
The barotropic primitive equation model on a global tropical 
mero 1S used to evaluates the pessibility of usamg @ numerical 
variational analysis as an initialization technique. 
Initialization of data with the nonlinear balance equation 


is used as a comparison. 
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i AN ERORUC oy 


Weather prediction in the tropics has always suffered, 
basically from a paucity of data. Hemispheric numerical 
merediction models not only suffer from this deficiency, but 
fomeine the boundary conditions in the tropics probably 
poo ceveriorates the product in this area. 

Jn recent years there has been some increase in 
Memeorological data over the tropical regions, particularly 
ia the form of commercial aircraft reports. This increase 
eiivowed FPiecet Numerical Weather Central, Monterey, to commence 
PomemjeCLive analysis for a global belt centered in the 
treaics. mm addition. FNWC was also investigating the use 
wee mumerical variational analysis scheme for this global 
belt. Since this scheme simultaneously solved for the 
wind and mass fields using the maximum amount of data, it 
could potentially be used as an initial state for a tropical 
fimerical prediction model, if the fields were dynamically 
consistent. 

Meebi-layer baroclinic numerical primitive equation 
Peeawetion models for limited regions in the tropics have 
been developed by Krishnamurti (1969), Miller (1969), 
Harrison (1969), and Elsberry and Harrison (1971). Vanderman 
and Collins (1967) have used a barotropic primitive equation 
model for operational experiments with data from a global 


meepical pelt. 





The availability of tropical data for 2 gilobal bella, 
allong with the limited-region prediction mode Gere loOpedmhy 
misberry and Harrison provided the impetus 10, scm Worl. 
mhe object of this thesis then, is to perform experiments 
involving the initialization of data and the development 
of a primitive equation model for short range tropical 


forecasts on an operational basis. 





II. PRIMITIVE EQUATION PREDICTION MODEL 


A. BASIC MODEL 

A ten-layer primitive equation model, developed by 
Harrison (1969) and Elsberry and Harrison (1970, 1971), 
Was used to formulate less complex models for possibile 
operational use. The original model was developed for short 
range predictions over a 34x15 point tropical grid, with a 
erid distance of 1.5 degrees longitude at 22.5 N. A 
mereui1onless, adiabatic version of this model required 
mapneximavely 20 minutes of computer time to produce a 
twenty-four hour prognosis. 
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The ten-layer model averaged parameters over two points 
only, in the direction of the space derivative. Finite 
differencing on a constant pressure surface, for the general 


variable S is as follows: 


Ne 


2, (sv) -{ | s(r, 144, K)| | v(I,S+1,K) aevenere® 


m(J+1) m(J ) 
-(s(z,5-3,x) | y ) 2 wsd=LE)]} 2 
m(J ) m(J-1) Ay 
where 
S(T,J+3,K) = SCLdtl.K) + S(1.d,K) | 
2 
ewan, J-3,K) = SOE a SEG) 
2 @ 
mn “meeebent foature cf the madel is the verticel distri- 


Duvron of variables, as shown in Fig. 1, since vertical 
Meelon 1S not averaged. All other parameters are averaged 
teerornly in the vertical, except at the lowest level, 


Vertical differencing is as follows: 


4. os) = { [w G.9.e+1)] | (1,5, K+1) | 


-[ (1,9,x-1)] [s(1,0,K-1)] ak 


where 


S(I,J,K+1) = S(1.5.K+2) + S(T,5,K) — 
2 
S(1,3,K-1) = S(2.3.K) + S(1.5.K-2) 
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Arakawa (1966) showed this finite difference flux form 
conserves the mean square of potential temperature and mean 
equare kinetic energy while preventing nondliimear instability 
caused by aliasing. The model remained computationally 
stable with minimal smoothing or averaging utilizing a six- 
peamute time step. 

The first time step is forward, where for the general 
parameter S, time step, n, and time difference, At, 

n=l = ee eso 
Centered or leapfrog time steps are used for subsequent 
Omme steps, where 

oe = Ss th At f(s”). 

Reinitialization, to avoid solution separation, during the 


-t. . ohn — Ts , - 3 i. « a ov =~ ~~ ere 
meee -oOur our Per1tog wee uo. Lewmd vo be necessary,. 


B. BAROTROPIC ONE-LAYER MODEL 

Modification of the ten-layer model into a one-layer 
barotropic model was the first step in the development of 
pm@eoperational model. The 250 mb surface was chosen as the 
Seamer Lorecast level due to the relatively high density 
@emedservacions in the tropics at this level, according to 
Krishnamurti (1970). 

Vertical distribution of variables in the one-layer 
Dbemeeuropic model is shown in Fig. 2. Required inputs are 
the winds and geopotential at 250 mb. A ficticious 1000 mb 
geopotential is derived, utilizing the barotropic assumption, 


— 


to maintain a similarity between the barotropic one-layer 


2 








and the ten-layer model. Since this) as @asbare tae nn. 
model Equations (3) and (4) are not required. In addition, 
the frictional terms in Equations (1) and (2) are neglected. 

Initial experiments were conducted using the 510 point 
grid that was used in the development of the ten-layer 
Mmoeel, This grid is illustrated in Fig. 5. Approximately 
two minutes of computer time is required for a twenty-four 
eur prognosis. 

Later work was performed using the FNWC (Monterey ) 
Tropical Global Mercator Projection, shown in Fig. 6. This 
fone 144x49 point grid with a grid distance of 2.5 degrees 
longitude at the equator. <A twenty-four hour barotropic 


prediction on the global band requires approximately 25 


Petes om.  oomniiter Game 
© ae wares ee ww . “y- : ae a 


C. BAROCLINIC TWO-LAYER MODEL 

mwa second step, a two-layer model Of the aenesphere 
meeeacveloped to include the effects of friction in the 
surface layer. With the exception of Equation (4), the 
Gull set of equations was used in this model. However, in 
BPaeteion to not predicting moisture, diabatic heating was 
neglected so that in Equation (3) only advection was per- 
Mrceoed to Change the initial temperature field. 

The vertical distribution of variables in the two-layer 
model is shown in Fig. 3. Required inputs are the winds, 
geopotential, and temperatures at 950 and 250 mb, and the 


geopotential and temperature at 1000 mb. 
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The friction terms are given vy the bulk aerodynamic 


transport relations 


a: 


BL, = ~g p) 
dp 
a 
a4 i dp 
where 
TT = g ca {V| u, 
1a == 3 cafV| vie 


mmeevalues of the drag coefficient, Cd, are 

Cd = 0.5 ly {2 x 10 79, 
for winds from 1 to 15 meters per second, and 

ea = 2.6 x 10°, 

TOr wind speeds greater than 15 meters per second. The 
od USSG in the above cvaiculation is a Jinear Timetion 

of the 950 mb wind with the direction held constant and 

the speed reduced 20 %. Friction is considered only in the 
Mager between 1000 and 900 mb, where the stress and the 
feeureal Sradient of stress are both assumed to be zero 

at 9OO mb. 

Pemsobaved Carlier, the upper level of most interest was 
eeewe50 mb level. However, this implies that the wind field 
mmmee very deep upper layer is defined by a level near the 
Mpper extreme of the layer. Consequently, the two layer 
baroclinic model was also tested in experiments with the 


350, 450, or 550 mb level replacing 250 mb as the upper 


Wind level. 


144 








All experiments with the two-layer model were conducted 
using the 510 point grid. A twenty-four hour prognosis 


required approximately seven minutes of computer time. 


D., BAROCLINIC THREE-LAYER MODEL 

Inclusion of data at the 750 mb level produced a three- 
layer version of the baroclinic model to provide a more 
realistic approximation of the atmosphere. The 800 mb 
above the surface layer was divided into two layers at the 
500 mb level. The development of this model also presented 
the opportunity to evaluate the effect of increased vertical 
feeoluvion versus the resulting increase in computer storage 
and time requirements. Fig. 4 illustrates the vertical 
distribution of variables for the baroclinic three-layer 
model. 

All experiments with this model were conducted with the 
Seeoepoint grid. Computer time of approximately nine minutes 


meomeeourred for a twenty-four hour prediction. 


EB. BOUNDARY CONDITIONS 
ine computationally stable integration of the primitive 
equations requires a set of boundary conditions that are 
Sememotent with the finite difference scheme. 
Flux through the top of the domain is eliminated by 
Pompine OG) = 0, 
At the lowest level, 1000 mb, Equation (5) is imposed. 
Experiments were conducted with three sets of horizontal 


boundary conditions. 








ia oat lusipoundanies 

Blsberry and Harrison (1971) discuss the boundary 
conditions which are applicable to the basic model. 
Initialization was accomplished using the nonlinear balance 
equation after computing the rotational component of the 
wind. 

The requirement for east-west boundary conditions 
on the 510 point grid is eliminated by making the data 
cyclical in the initialization. Naturally, east-west 


boundary conditions are not necessary on the global band. 


By setting 
m+ 
ee - “oe ye, Ws. 5 = 2uU,, aod 14 = 
Jagik ae ea aad Jib Psa ? Shes 
mJ 
feeee oo. , 2s Che norunern=mest row in the gvidyg no flax of 
ed. - 


mass or energy is allowed. The southern boundary is treated 
aie Similar manner. 
ee Restoration or constant-flux boundaries 

Kesel and Winninghoff (1970) describe a constant- 
fiux, restoration boundary technique that eliminates both 
the necessity for altering the initial mass structure and 
the faise reflection of the computational mode at the 
boundary. Initial height and wind fields near the boundaries 
are preserved. At the completion of each time step the 
newly predicted values of the fields are restored toward 
values at the previous time step with a specified field of 
restoration coefficients. Figure (7) illustrates the 


restoration field near the northern boundary. The procedure 


16 








results in a fully dynamic forecast seven rows in from the 
Boundary, a persistence forecast for the two rows at the 
boundary, and a blend in-between. Restoration boundary 
Memeitions were restricted to experiments using the global 
meid. 
Be. Interpolated boundary conditions 

Krishnamurti (1969) used implicit north and south 
poles to define boundary values. Geopotential values at 
M@eenorth pole are defined by the mean value at row Jj, 
Where the north boundary is row j+l. The three components 
Oi the wind are assumed to be zero at the pole. A linear 
Mterpolation between the pole and row j is performed each 
mame Step to calculate the values of the four parameters 
me row jit. A Simi specifics the southern 
moumcdary. Kkxperiments with interpolated boundary conditions 


were confined to the global band. 


di 








Tit. DATA ANITEALIZATION 


hae importance of initial conditions for a primitive 
equation model has long been recognized. Raw meteorological 
Maca, used as an input, guickly results in high frequency 
inertio-gravity oscillations. Shuman and Hovermale (1969), 
Benwell and Bretherton (1968), Smagorinsky et al. (1965), 
Charney (1955), and Hinkelmann (1951) all propose an initial 
eos With no divergence or vertical motion, to suppress 
these oscillations. Phillips (1960), however, showed that 
Bor 2 Simple linear model an initial divergence equal to 
Maat SUgeested by quasi-geostrophic balance is required. 
Hinkelmann (1959) demonstrated that high frequency 
SBscillations will be filtered out if the time derivative 
emeeerereence vanishes in the initial field, that is if 

z2.- 0 (9) 


at” 
Meee 1 is the divergence and n=1,2,3,.....- Note that 





Meine the balance equation to specify the initial conditions 
Mesuits in satisfying the above conditions for n=l. Finally, 
this condition has been satisfied for n=1 and 2 by Miyakoda 
and Moyer (1968) and Nitta and Hovermale (1969), by 
integrating the primitive equations forward and backward 
Micil the initial field meets the condition. Two types of 


initialization were tested in this series of experiments. 
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A, NONLINEAR BALANCE EQUATION 

Elsberry (1969) utilized the nonlinear balance equation 
and special boundary conditions to formulate a diagnostic 
model that would match the prognostic model (ten-layer with 
no-flux boundary conditions). The nondivergent part of the 
wind was obtained from 

V “ys 7, 

where 1d is the relative vorticity from the observed wind 
memoonents u and v. Smoothing the relative vorticity fields 
mace vertical improved the vertical consistency of the 
stream-function fields. No normal component of the wind was 
allowed at the north and south boundaries, implying that 
these boundaries were streamlines. Thus ¥n and Y s were 


m oy oe ~ As -t. ~ as BL +. _ 8 tA in 4 i 
Sec equal to constenbes 2t each af the bounsaries with 


a . + =-udy, 


where u is the mean zonal wind component averaged over the 
erid, and Ay is the distance between the north and south 
boundaries. In the manner of Krishnamurti (1969) the grid 
was extended by five points between the east and west ends 
meeecone Srid which are jinearly matched, providing cyclical 
Pomuinuity. 

The Yfield was used to eee the geopotential from 


mame Nonlinear balance equation, 


24=V- oY oY 
Vega Virvyser (gh ot) 


Boundary conditions on ® at the northern and southern walls 


were made geostrophically consistent with the streamfunction 


ALY. 








cen Lerence aeaace by using a mean value of f across the 
limited region grid. On the global belt the averages of 
the observed geopotential heights along the walls were used for 


the boundary conditions. 


B. NUMERICAL VARIATIONAL ANALYSIS 

Sasaki (1958) introduced the variational method to 
Spcain dynamically consistent initial fields, showing 
theoretically that quasi-geostrophically balanced or 
Givergence-free fields could be obtained as a solution to 
See Huler-Lagrange equations. Filtering high frequency 
oscillations using Equation (9) as a low pass filter, where 
Weis nov limited to divergence, but may be any meteorological 
parameter to be filtered, was also proposed by Sasaki (1970). 
hong wave equavions are usea with quasi-svecary 
ditions as a dynamical constraint in a numerical variational 
analysis. This method, as adapted by Lewis (1971) for 
operational use by FNWC, Monterey, was tested as an 
Meet tiZzaction scheme for the barotropic prediction experi- 
Memcs On the giobal grid. 

The variational analysis was obtained by minimizing the 


functional 
5 = ff {2 [(u-2) + (v-¥)?] + F (2-7)* 
S 
+ oc (28)° ee (2)? } as 
where the integration is over the surface S being analyzed. 
Here U, v, end z are the observations interpolated to grid- 
points by an objective analysis scheme, and u, v, and z 


are the desired fields which result from the minimization 


20) 








of J, and ey are weighting functions dictated by the 
density and reliability of observations. The last two terms 
in the functional can be interpreted as a control on the 
amount of acceleration in the analysis, and are related to 
the low pass filter characteristic mentioned earlier. The 
larger the ratio of the dynamical weight % to the obser- 
VabLonal weights, e and ca the more the analysis will be 


forced toward a steady state where 


O 


pu Cll OZ ly. 
me - 5 (Dx Y-Vu + fv 


and 


My. 1 92 wiv - Ff 
3 Sy Y-Vv UW 
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It 


Q. 


ct 


pmeady State can only be attained in the limit, if oY, 


— 
ee 


/ 
. = sh ms = nts TwWe 4 Vrs TAs a; 
weer yy are Fer Gc Tee ol ae ae 


OUOv ei SReMWeIe frend. ; 
peme Small tendencies in the fields however, may be desirable 
mere tne initialization of a primitive equation prediction 
mode | . 

Ho enhance the operational value of the variational 
analysis scheme, Lewis (1971) developed a technique to 
avoid solving the three equations for the u, v, and z fields 
Simultaneously. However, this technique involved grouping 
the inertial terms with the local change terms, and minimizing 
the imbalance between the remaining terms, the pressure 
aeadi1ent and coriolis terms. The usefulness of the variational 
analysis for initialization was degraded by this technique, 
Pervicularly in high wind regions. The technique was then 
modified to include an approximation of the inertial terms, 


by evaluating them using the input fields, 0 and V¥. 


call 
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A, PREDICTION EXPERIMENTS WITH THE LIMITED GRID 

Mes group of experiments is an attempt to test how much 
a barotropic model and simple baroclinic models degrade the 
results of the ten-layer model. Data initialization was 
meecomplished using the nonlinear balance equation technique 
and the prediction models had no-flux boundary conditions. 
Mie input data contained both an upper level cold low and 
mee anticyclonic flow over an incipient tropical storm. 

ime Barotropic Model 

The barotropic model provided a twenty-four hour 
mpemees very favorably with the tweamy-four 
hour results from the ten-layer model. Figures 8 and 9 
miiustrate the twenty-four hour prognoses from the respec- 
tive models. The barotropic prognosis shows essentially 
mae same height pattern as the ten-layer prognosis. The 
M@eeuest difference is in the height of the cold low, 

Meere the barotropic prediction is approximately ten meters 
higher than the ten-layer prediction. 

Figure 10 illustrates the change of the total, 
monal, and eddy kinetic See ee With respect to time fer 
the barotropic model. Although an apparently spurious 
@eeecase in the eddy kinetic energy is noted, no degradation 


of the twenty-four hour prognosis was noted. 


Ge 








Cee uM@elayer Model 
The attempt to predict the 250 mb surface with the 


mro-layer model was not successful. The exponential 
srowth of eddy kinetic energy associated with this failure 
meeshnown in Fig. 11, where the solid lines indicate the 
250 mb kinetic energies. The dashed lines illustrate the 
450 mb kinetic energies. Two-layer models with upper levels 
meee s5O0 and 550 mb also did not exhibit unreasonable growth 
Meeecne eddy kinetic energy. This implies that perhaps 
euecessful predictions could be carried out when the broad 
moper layer is defined with data from levels lower than 
250 mb, or with input data which has been averaged over 
meveral layers. 
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Twenty-four hour prognoses from the three-layer 
meecl were very encouraging. The twenty-four hour 250 mb 
meoenosis, as shown in Fig. 12, is almost identical to the 
memetayer prognosis shown in Pig. 9. There are no appreciable 
mmerercnces in the heights, Of course, the three-layer 
Meeel also provides predictions at 750, 950, and 1000 mb. 
The (50 mb twenty-four hour prognoses as illustrated in 
Pag. 13 and 14 also have no significant differences. Fig. 15, 
the three-layer 1000 mb Creer oue NOU OeOemMOcds.- Wiltce 
having the same pattern as Fig. 16, the ten-layer version, 
@eeetilied the tropical cyclone approximately ten meters 


Mere, thus reducing the gradient surrounding the storm. 


ae 








Thus a three-layer model which would incorporate 
both the boundary layer and the upper tropospheric flow, 
which are known to be important in the tropics, may be 
measible. The middle layer would be a useful buffer and 
mao an additional degree of freedom for baroclinic develocp- 
Ment. tthe scarcity of data in the tropics for this mid- 
mayer could lead to the use of low cloud motions, as 
measured from the Applied Technology Satellite, for 
mieoialization of this level. 

Unfortunately, only one data set was available for 
mae experiments with the limited grid and further experi- 
ments with more data sets are required before these results 


may De considered conclusive. 


fee EXPERIMENTS ON Tess GLOBAL BAND 

Use of the nonlinear balance equation for the 
Merotalization of the primitive equation model has been 
Bemonstrated previously by Harrison (1969) and Eilsberry and 
Harrison (1971). Krishnamurti (1969) also used this tech- 
meeme Lor the initialization of a primitive equation model. 
Bae predictions obtained from the barotropic and two- and 
three-layer models tend to confirm these results. Although 
only the rotational component of the wind and a smooth height 
field from the balance equation are used as initial data, 
@eeealistic prediction of systems is the resultant. 

It was considered sufficient to test the initialization 


properties of the numerical variational analysis with the 


eu 








merotropic model in order to save computer time. To assure 
m@ cood data set, climatological fields were used rather 

than actual data in the global band experiments. Both 500 mb 
me 500 mb climatological fields were used. Although only 
mie 500 mb results are discussed and illustrated, the same 
mepe problems were encountered with both fields and the 
Memacions were similar. Since the size of the FNWC tropical 
Mercator projection would not permit a detailed view of the 
entire grid in this publication, only a small portion of 

Sie serid is presented. This portion was chosen as being 
Bepresentative of the problems and solutions encountered 
With the whole grid and includes the entire north-south 
Sxtent of the grid, the western Pacific Ocean, the Asian 
Beastiine, Japan, Indonesia and Australia. The 500 mb 
meimacological field for January in this area features a 
broad trough in the north, a band of maximum winds between 
SO and 40N, flat high pressure from 20N to 20S and basically 
zonal flow in the southern hemisphere. 

Pigure 1/f represents the nonlinear balance equation 
Pbolalization, Solid lines are the initial data while the 
meemea Jines are the departure of the initial data from 
the 500 mb climatology. Smoothing near the boundaries is 
meeerent as is a large overall height rise caused by the 
Pepcialization boundary conditions. On a limited tropical 
mena, setting the heights on the north and south boundaries 
meme! to the mean height is reasonable, however, on the 


global band, with boundaries well into mid-latitudes, this 








meno Longer true. In the region of a trough allong the 
poundary, the field will be relaxed toward a boundary that 
mmmroo Nish, while the opposite effect will occur in the 
area where a ridge is along the boundary. New boundary 
conditions will be required if the nonlinear balance 
Sewaoion initialization technique is to be used on a grid 
mepe boundaries in mid-latitudes. 

Initial heights as determined by the numerical variational 
analysis technique are shown in Fig. 18 as solid lines. The 
Meenead Jines depict the deviation of the initial heights 
imeem the climatological 500 mb data. It is readily apparent 
@aeo this technique does not significantly alter the input 
iPreld. 

rernaps the greatest advantage to be realized by using 
eae numerical variational analysis technique for initialization 
a8 provided by the opportunity to vary the weighting functions 
Syer the analysis. itn sparse data regicns the dynamical 
Meemetiy, & , can be increased. In mid-latitudes the height 
meld weight, B ,) Cole same reasicds swallow ii one = uaeimles 
the wind field weight, > cam be -inereased. ~The proper 
Utilization of the weighting functions should guarantee a 
dynamically consistent input for primitive equation models. 

Lo determine if the variational analysis technique 
would initialize the barotropic model, a twenty-four hour 
Prediction with initialization by this scheme was compared 
tO a Similar prediction with initialization by the nonlinear 


balance equation technique. Restoration boundary conditions 
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were used in the model with the variational initialization. 
Figure 19 illustrates the twenty-four hour prediction for 
mis technique. The twenty-four hour prediction for the 
model initialized with the balance equation is shown in 
mee. 20. This model used no-flux boundary conditions. 
Menmoring the height differences caused by the initialization 
techniques, and particularly in regions away from the 
meundaeries, the predictions are very similar. Oniy the 
developing short wave features in the wind maximum between 
30N and 40ON are disturbing. In ali, three sets of boundary 
Conditions, no-flux, restoration, and interpolated, were 
tested with the two initialization techniques. All combi- 


nations developed the short wave features mentioned above. 


Rha nredanmnin 
we fl ee eet 


am 


antiy zanal southern 
m_mesomere, nor were there any short period height oscillations 
at the equator as noted by Winninghoff (1971) and Krish- 
namurti (personal communication). The conclusion to be 

Beewn from these limited experiments then is that the short 
Seemewieatures are independent of either boundary conditions 

Or initialization and that the variational analysis technique 
mer cuccessfully initialize the barotropic primitive 

Beuwev2on prediction model. 

The solid curves in Fig. 21 illustrate kinetic energy 
Beenees for the restoration-variational analysis prediction 
mece! for the forecast period. The growth of kinetic energy 
is attributed to an accumulation of energy in the smallest 


meomitted scales as the eddy kinetic energy appears to be 
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mie torcing function. This phenomena was also noted by 
Kesel (1968), in a barotropic primitive equation experiment 
mes a hemispheric grid. 

The accumulation of the eddy Kinetic energy was controlled 
by adding a term to the momentum equations to simulate 


meperal diffusion of momentum, in the form 
2 a 
ye m Vou and )Cm Wo V5 


where )(, the tavera leddy= cell us Vem ecm n Caer toed. 
meme oon of the grid distance. 

A twenty-four hour prognosis as obtained from the 
meeouropic model, with lateral diffusion terms included 
mt=5 x 10°), restoration boundary conditions and 
Bereacional analysis initialization is shown in Fig. 2e. 
Neve that the aforementioned short-wave features are absent. 
Seemmeconed lines in Fig. 21 illustrate the Kinetic energy 
meemeses for the model with the diffusion term included. 

Miat the lateral diffusion term prevents accumulation of 


meeeey in the small grid scales is evident. 
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Ve. CONCLUSIONS AND RECOMMENDATIONS 


pe barotropic or three-layer primitive equavion vredicpien 
model for use in the tropics on an operational basis seems 
well within reach, whether on a limited region or a global 
basis. The basic model for use on the FNWC tropical 
Mercator grid has been tested. Further experiments are 
required to isolate the cause of and either eliminate or 
Control the accumulation of eddy Kinetic energy. As the 
prediction models become more refined, a staggered grid 
System should be investigated to further reduce computer 
time and storage requirements. 

The suecessful use of the nonlinear balance equation for 
initialization of a limited-region primitive equation 
prediction model has been demonstrated by Krishnamurti (1969) 
and Elsberry (1969). However, if this scheme is to be 
utilized when the boundaries are in areas with height gra- 
dients larger than those encountered in the tropics, more 
BulLcable boundary conditions must be developed. Over a 
global band, no-flux boundary conditions may not be crucial 
and other boundary conditions may be applicable. 

The numerical variational analysis technique has been 
meepveed for use in initializing the barotropic primitive 
equation model. Using the nonlinear balance equation 
Initialization scheme as a comparison, it would appear as 


if the numerical variational analysis will successfully 
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in the two-layer model. 
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Fig. 4. Vertical distribution of variables 
in the three-layer model. 
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Pig. {. Restoration coefficients near the 
northern Doundary of the tropical global 
belt. A persistence forecast results at 
the boundary, a fully dynamic forecast 
seven rows in, and a blend in-between. 
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Total (T), zonal (Z), and eddy (E) kinetic 
energies vs. time for the barotropic model 
on the limited grid. 
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Fig. 11. Total (T), zonal (Z), and eddy (E) 
Kinetic energies vs. time for the two- 
layer model on the limited grid. Solid 
lines indicate 250 mb data while dashed 
lines indicate 450 mb data. 
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Fig. 17. Initialization of the 500 mb height field by the 
nonlinear balance equation technique. Solid 
lines are initial data while dashed lines are 
the departure from the 500 mb climatology used 
as input. 
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Pipe. 


val aaoLOnwyOt ENnewaeO My Neier telaey, te 
numerical variational analysis technique. Solid 
lines are initial data while dashed lines are the 
departure from the 500 mb climatology used as input. 
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Twenty-four hour 250 mb height forecast from the 
barotropic model on the global grid with restoration 
boundary conditions and numerical variational 
analysis initialization. 
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Fig. 20. Twenty-four hour 250 mb height forecast from the 
barotropic model on the global grid with no-flux 
boundary conditions and nonlinear balance 


equation initialization. 


one 








100 


Us 





O 6 iz 18 C4 
FOURS 


ie. 21. Total (T), zonal (Z), and eddy (E) kinetic 
energies vs. time for the barotropic model on the 
global grid. Solid lines are for the basic model 
while dashed lines are for the model with diffusion 
terms included. 
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Twenty-four hour 250 mb height forecast from the 
barotropic model with diffusion terms, restoration 
boundary conditions and numerical variational 
enigasets dy ra lizaclion. 
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